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Abstract 

We argue that the relative thermal conductance between interfaces with different morphologies is 
controlled by crystal structure through M m i n /M c > 1, the ratio between the minimum mode count 
on either side M m , n , and the conserving modes M c that preserve phonon momentum transverse to 
the interface. Junctions with an added homogenous layer, “uniform”, and “abrupt” junctions are 
limited to M c while junctions with interfacial disorder, “mixed”, exploit the expansion of mode 
spectrum to M m i n . In our studies with cubic crystals, the largest enhancement of conductance 
from “abrupt” to “mixed” interfaces seems to be correlated with the emergence of voids in the 
conserving modes, where M c = 0. Such voids typically arise when the interlayer coupling is weakly 
dispersive, making the bands shift rigidly with momentum. Interfacial mixing also increases alloy 
scattering, which reduces conductance in opposition with the mode spectrum expansion. Thus the 
conductance across a “mixed” junction does not always increase relative to that at a “uniform” 
interface. 
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FIG. 1. a) Abrupt interface, b) Interface with random atomic mixing (mixed interface), c) Interface 
with an added homogeneous layer at the junction (uniform interface). 

I. INTRODUCTION 


For over half a century, the thermal energy flow across solid-solid interfaces has been 
studied with only partial understanding of the underlying processes jl-3]. A microscopic un¬ 
derstanding of these interfacial thermal processes requires deconstructing thermal interfacial 
conductance, which brings many challenges, including consideration of a broad spectrum of 
interacting dispersive phonons, varying mean free paths, and additional phonon interactions 
with defects, impurities and other interfacial imperfections jd]. Moreover, as the spacing be¬ 
tween two interfaces reduces to distances on the order of the phonon coherence length, wave 
interference and coherent transport contribute to the thermal resistance in a non-additive 


fashion 


8J. 


Early models of interfacial thermal conductance focused on the effect of acoustic matching 
9], nonlinear dispersion [loj], and bonding 101 on perfectly abrupt interfaces (Fig. QJi). 
Interfacial imperfections were later included in the diffuse mismatch model (DMM) as sources 
of diffuse scattering jll]. Although this model is widely used, it does not accoimt for 


atomistic interfacial details 12], which have been shown to affect interface conductance 
measurements 4 ] • 

One interfacial imperfection is random atomic mixing (Fig. [Dd), which can be a frequent 
byproduct of nanostructure fabrication. The addition of random atoms at an abrupt in- 
ace generates two effects in the harmonic regime: 1) it changes phonon transmission 
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13dl6|: and, 2) it couples phonons with different transverse wavevectors ( k ±) by breaking 


the translational symmetry at the interfacial plane 


17]. Some papers focused on the effect 


of mixing on transmission (Effect 1) and showed the importance of the frequency depen- 
, the correlation length of the random distribution [mI, 15] and the acoustic-optic 
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coupling 


16]. 


Interestingly, the papers that focused on the coupling of phonons (Effect 2) noted that 
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the thermal conductance of the mixed interface was larger than that of the abrupt interface 
for a simple cubic crystal interface [l7, 18] and for a Si/heavy-Si interface 19]. Moreover, 


Kechrakos 


181 ] noticed that the conductance of the mixed interface was even larger than the 


conductance of an interface with an added homogeneous atomic layer (Fig.QJ), which we call 
a uniform interface. Those results suggest that adding disorder at interfaces increases con¬ 
ductance, contrary to bulk materials where adding disorder has the opposite effect. However, 
the role of crystal structure on the conductance of the abrupt, mixed and uniform interfaces 
remains unclear. 

In this paper, we demonstrate that the relative conductance between interfaces with 
different morphologies (Fig. [T]) depends on crystal structure. For the systems considered in 
this study, adding a unit cell monolayer of mixing to an abrupt interface always enhances 
the interfacial conductance, but the extent varies over an order of magnitude according to 
the crystal structure (Fig. [2]). Furthermore, the conductance across a mixed interface does 
not always increase relative to that at a uniform interface. In fact, while the conductance 
increases for simple cubic (SC) and diamond cubic (DC) crystal structures, it decreases 
for face centered cubic (FCC) crystals. This suggests that the commonly invoked virtual 
crystal approximation, which models the mixed interface as a uniform interface, alternatively 
overestimates or underestimates the effect of interfacial mixing on thermal conductance. For 
DC crystalline interfacial regions, we show that the enhancement of conductance by mixing 
depends on phonon polarization. For instance, mixing increases transmission between TA¬ 
TA modes but not between LA-LA modes. 

We explain our results within the Landauer theory, where thermal conductance is directly 
related to the product of the number of conducting channels or modes (M) times their 
average transmission (T). We find that 1) the crystal structure determines the relative 
magnitude of the minimum of the contacts’ modes M min vs. the conserving modes M c that 
conserve the component of phonon momentum transverse to the interface. On the other 
hand, 2) the interface morphology determines if phonons move through M m i n for mixed 
interfaces, or M c otherwise. Based on these two concepts, we show that the conductance 
across a mixed interface increases relative to that at a uniform interface when M m i n > M c , 
with larger degree of enhancement as the inequality increases. The larger enhancements, 
seen in SC and in TA branches of DC, are associated with the emergence of voids in the 
conserving modes (M c = 0). Such voids in turn arise when the subbands shift but do not 
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distort with increasing momenta. 

We start by deriving an inequality (Eq. [5]) between the conductance of the mixed and 
uniform interfaces from the Landauer theory (Sec. [HI) - Then, we describe how the modes 
(Sec. IIIII) and transmission fSec. UVD shape that inequality according to phonon polarization. 
In Sec. IIV1 we derive analytical expressions for the transmissions of the scalar SC and FCC 
systems. For the uniform interface we find a maximum conductance when the junction mass 
is the arithmetic mean of the contact masses. For the mixed interface, we find that the 
transmission between phonons that do not conserve transverse wavevector, k±, depends on 
the difference of the contact masses and on the alloy scattering factor, a(l — a) with a the 
fraction of heavy atoms at the interface. 


II. LANDAUER DESCRIPTION 


Thermal conductance G q is defined as the ratio between heat flux q and temperature 
drop AT. Within the Landauer theory this quantity can be expressed as 20] 


G q = 


q 


i q 


AT AAT 


oo 

1 fdco^ dN , _ 
— / — fuj——MT , 
A J 2tt dT 

o 


( 1 ) 


where I q is the heat current, A is the cross-sectional area, hu is the energy carried by a 
phonon, N is the Bose-Einstein distribution, M is the number of propagating modes, which 
we refer as “modes” throughout this paper, and T is the average transmission per mode. For 
a given contact and frequency ui, the propagating modes are the eigenvectors (x n oc e d k ^n-^t)^ 
of the equation of motion for the contact with eigenvalue cu 2 , with real wavevector k and 
with group velocity in the transport direction. The product MT equals the sum of the 
phonon transmissions between modes on the left and right contacts. This quantity can be 
calculated from non-equilibrium Green’s functions (NEGF) as MT = Trace{T;GT r Gd}, with 
G the retarded Green’s function and T the broadening matrix for the left (/) and right (r) 


contacts 


21 


23]. 


For the uniform interface (Fig.QJ), the symmetry in the transverse direction requires that 
phonons crossing it conserve their transverse wavevector k±. Thus, the nonzero contributions 
to MT are transmissions T k , k , between contact modes with the same k±. Referring to the 
number of these transmissions as M c , the conserving modes, and their average as T c , we can 
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express MT for the uniform interface as 


MT uni = J2T k± , k± = M C T C < M c . 


( 2 ) 


M c is given by the overlap between the projections of the frequency isosurfaces of the contacts 
onto the k± plane (Fig [3b). Note its role as an upper bound of MT uni . Also note that the 
abrupt interface (Fig. [TJl) is a limiting case of the uniform interface. 

When we replace the homogeneous interfacial layer of the uniform interface by random 
contact atoms (Fig. [TJd), the atomic disorder breaks the transverse symmetry and allows 
phonon transmission T k±jk > ± between modes that do not conserve k± 17]. That disorder also 
decreases the transmission between modes that conserve k±. We can express MT for 

the mixed interface as 


MT mix — T k±tk 

ki k±^k'. 


— M C (T C — 5T C ±) + M nc 5T nc ^ < M n 

non-conserving 


v 

conserving 


( 3 ) 

( 4 ) 


where M nc 5T nc -|. represents the increase in conductance due to the newly available channels. 
Note that energy conservation ensures that MT mix is bounded by the minimum of the bulk 
modes of the contacts: M min = min(Mg M r ). Thus we define M nc = M min . 

Comparing Eq. [2] and [4j the conductance of the mixed interface is larger than that of the 
uniform interface, G m i X > G un i , if 



duj n dN M AT 1 
—hu)——M m i n dl nc -[• 

27t a T 


> 



2vr dT ^ 


( 5 ) 


In other words, G mix > G uni if the gain in conductance by opening new channels that do 
not conserve k±_ (M nc 5T nc - 1 -) surpasses the loss in conductance by phonons conserving k± 
(. M C 5T C ±) over a window set by the cut-off frequency and the temperature. 


III. MINIMUM VS. CONSERVING MODES 

We calculate the harmonic conductance of abrupt, uniform and mixed interfaces embed¬ 
ded into four different crystal structures: 1) SC and 2) FCC crystals, where the atomic 
movements are simplified to a single direction and thus the interatomic force constants 
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(IFCs) are scalars; 3) FCC crystal with IFCs calculated from the Lennard-Jones (LJ) po¬ 
tential; and 4) DC crystal with IFCs calculated from density functional theory (DFT). The 
interfacial region for each system consists of a monolayer of primitive unit cells (Fig. 1). 
The same IFCs and lattice constants are used throughout each system to isolate the effect of 
mass disorder, which has been proven to dominate the scattering of cross-species interactions 


24]. The ratio between the atomic masses of the contacts is 3 for the SC and FCC systems 


and 2.6 for the DC system, corresponding to the mass ratio of Si and Ge. The conductance 
is calculated using NEGF, and the details of the simulations and assumptions are given in 
Appendix [A] 



FIG. 2. Plot of the average MT normalized by the results for an abrupt interface versus the average 
mass at the junction layer for the cases described at the beginning of the section. Adding mixing 
at an abrupt interface enhances the interfacial conductance in all systems simulated in this work, 
but the extent depends on crystal structure. Compared to uniform interfaces, however, mixing 
does not always yield an increase in conductance. 

Figure [2] plots the frequency average of MT, which one can interpret as conductance 
without the “low” temperature dependence (Eq. [T]), vs. the average mass at the interfacial 
layer (mj) for the different crystal structures. For these systems, adding a unit cell monolayer 
of mixing to the abrupt interface always enhances the interfacial conductance, but the extent 
depends on crystal structure. For the scalar SC crystal (Fig. [2jr), the maximum relative 
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change of ( MT) between the mixed and abrupt interfaces is about ((2.5 — 1)/1 x 100%) = 
150%. However, for the FCC and DC crystals (Fig. |2]o, c and d), the relative change of 
(MT) is only about 13%. Furthermore, the conductance across a mixed interface does not 
always increase relative to that at a uniform interface. In fact, it increases for SC and DC 
crystals but decreases for FCC crystals. 

The increment of conductance from the abrupt interface to the mixed interface (Fig. [2]) 
relies on the atomic extent of the mixing region. For this special case, Eq. ED tells us that the 
gain in conductance by phonons that do not conserve k± surpasses the loss in conductance by 
phonons conserving k± ( M nc ST nc ^ > M C ST C jJ. As the extent of the mixing region increases, 
phonon back scattering increases and transmission decreases. Thus, 5T nc ^ decreases while 
5T C ± increases, making the inequality more difficult to be satisfied. At some point, the 
inequality stop being true and G m i X becomes less than G a b r , which is the usual experimental 
outcome |4|. 



W=2.7xl0 13 w=3.1xl0 13 cj=4.2x10 13 



FIG. 3. a) For SC, the conductance for the mixed interface is larger than the uniform interface 
because the loss in transmission, T, due to the extra scattering brought by the random atoms (area 
with vertical lines) is dominated by the gain in modes, M, coming from transmissions between 
phonon that do not conserve k± (area with horizontal lines), b) While the MT un i spectrum 
is limited by transmissions between phonon conserving k± (overlap region), the extra MT m i x 
spectrum comes from transmissions between phonons that do not conserve k±. 

For the SC crystal, the large conductance increase of the mixed interface results from the 
wider MT spectrum (Fig. [3ja) . This extra spectrum comes only from transmissions between 
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modes that do not conserve k±. In fact, over that frequency interval, the available contact 
modes do not overlap (Fig. [3)o) , banning transmissions between modes conserving k±. Thus, 
M c = 0 and MT uni = 0. Mixing removes the requirement of conserving k±, opening M min 
conduction channels and making MT mix > 0. 

Figure Ek shows MT for the mixed, uniform and abrupt interfaces. MT mix is split into 
the contributions from modes conserving k± and those that do not. This gives us a pictorial 
representation of Eq. [5j G m i X > G un i because the MT area gained due to transmissions 
between modes not conserving k± is larger than the MT area lost due to disorder among 
the modes that conserve kj_. 



b) Available modes in the contacts 

oj=4.5x10 13 w=5.5x!0 13 w=6.0xl0 13 



k x a k x a k x a 

\M min =M c \ \M min =M c \ \M,ni„>M c \ 


FIG. 4. a) For FCC, the conductance for the uniform interface is larger than the mixed interface 
because the loss in T due to the extra scattering brought by the random atoms (area with vertical 
lines) dominates the gain in M coming from transmissions between phonon that do not conserve 
perpendicular momentum (area with horizontal lines), b) M c ~ M m i n over most of the spectrum. 
Thus, mixing provides little advantage by allowing transmission between modes not conserving k±. 

A similar pictorial representation for the scalar FCC interfaces is shown in Fig. [4^. In 
this case, G un% > G rmx because the MT area gained due to transmissions between modes 
not conserving kj_ is less than the MT area lost due to disorder among the modes that 
conserve k±. Note that MT uni and MT mix cover the same frequency range, and the overlap 
of the contacts’ modes M c equals their minimum M min over most of the spectrum (Fig. (Do). 
Thus, the accessible modes on the mixed interface M min do not bring any advantage over 
the existing modes M c on the uniform interface (Eq. [5]). The dominant conductance is then 



























































decided by the transmission, which in this case favors the loss in the conserving modes over 
the gain in the non-conserving ones. 

From Fig. [3] and SJ we note that the relative magnitude between M min and M c plays 
an important role determining the larger MT between the mixed and uniform interfaces 
(Eq. EJ). This is not surprising because of their roles as MT upper bounds for the mixed 
and uniform cases respectively. We can distinguish three cases: when 1) M c & M min , 
the modes conserving k± reach the physical limit of modes that can carry heat in one of 
the contacts. Equation 0 tells us that the transmission alone decides the dominant MT, 
which can be either the uniform or mixed MT. For the scalar SC and FCC structures, 
whenever M c « M min we see that MT uni > MT mix (Fig. [3] and SJ). Therefore the loss 
in transmission on the conserving modes surpasses the gain in transmission on the non¬ 
conserving modes (5T C ± > 5T ne j- in Eq. [5]). When 2) M min > M c , the dominant MT results 
from a balance between the added modes that do not conserve k± and the loss in transmission 
on the modes that conserve k± (Eq. E3). For instance, in the SC structure, MT mix becomes 
larger than MT uni as the ratio M min /M c increases (Fig. E}. When 3) M min > M c = 0, 
MT mix > MT uni = 0 and MT mix is only due to transmissions between modes that do not 
conserve k± as shown by Fig. El These three criteria may help in the search for interfacial 
materials where a particular outcome is expected from atomic mixing over the harmonic 
regime. 

For the diamond crystal, the polarization of the incident and transmitted phonons plays 
an important role in deciding the outcome of the dominant conductance. We give a brief 
description of the calculation process in Appendix |B} Figure [5^ shows that G mix > G uni 
mostly because of TA phonons in the light contact transitioning to TA and LA phonons in the 
heavy material. The transmissions between phonons that do not change polarization behave 
similarly to the scalar crystals. For TA-TA and LA-LA transmissions, MT mix > MT uni when 
M min > M c and MT uni > MT mix when M min M c (Fig. [5h>). The transmission between 
other polarizations will be analyzed in the next section. 

The polarized modes for the diamond crystal uncover an interesting similarity between 
the modes of the SC and the TA branches in DC and between the modes of the scalar 
FCC and the LA branch in DC (Fig. ED- For the cases where mixing significantly enhances 
conductance (SC and TA-TA branches in DC), we see a common central void in the modes 
for the heavy contact. This void arises in phonon bands where the k± ~ 0 subbands only 
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a) MT vs u per polarization 



FIG. 5. a) MT for DC crystal split in polarizations. The conductance of the mixed interface 
dominates mostly because mixing favors transmissions between TA-TA and TA-LA modes, where 
the ratio of M m i n /M c is larger, b) M m i n and M c for the different polarization branches. 

cover a fraction of the whole band spectrum (Fig. [7]). Indeed, after the cutoff frequency of 
those subbands, the k± & 0 modes, or central modes, start to become unavailable. From 
another point of view, the void originates when the upward shift of the kj_ ~ 0 subbands 
dominate their shrinking as \k±\ increases. We see this happening for SC and TA-TA but 
not for FCC and LA-LA (Fig. [7]). 

For the scalar SC and FCC crystals, the existence of the void can be associated with the 
independence of interlayer coupling as k± increases. For a SC crystal with atomic mass m 
and interatomic force constant /, the subbands are given by 

w 2 m = fon - 2faff cos(fc 2 a), (6) 

with the onsite coupling f on = 6/ — 2/ cos(k x a) — 2/ cos(k y a) representing the atomic inter¬ 
actions within a transverse layer of atoms, and the offsite coupling f Q ff = f representing the 
interaction between layers. As the magnitude of k± increases f on increases, shifting upwards 
the subband but f a ff remains constant keeping their width stable. On the other hand, for 
a FCC crystal the subbands are given by 

uj 2 m = fon -2 f off cos (k z ^j , (7) 

with f on = 12f — 4f[cos(k x a/2)cos(k y a/2)] and f 0 ff = 2f[cos(k x a/2 ) + cos(k y a/2)]. As 
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b) Available LA modes for Si and Ge 

oj=3.0x10 13 m=3.Txl0 13 w=4.0xl0 13 u=4.3x!0 13 



FIG. 6. Available modes in the DC contacts for the (a) TA and (b) LA branches. The modes for 
SC crystal and the TA branches, where mixing enhances the most the conductance of an abrupt 
interface, present a central void that enhances the ratio M m i n /M c . 


the magnitude of k± increases, f on increases shifting the subbands upward, but at the same 
time, foff decreases shrinking their bandwidth. 


scalar SC scalar FCC 



Mixing significantly Mixing marginally 
enhances G enhances G 



Mixing significantly Mixing marginally 
enhances G enhances G 


FIG. 7. When the shift in the subbands dominate their distortion, the k± = 0 subband only covers 
a fraction of the whole band spectrum. This generates a void in the conserving modes, which is 
seen when mixing significantly enhances the conductance of an abrupt interface. The solid lines 
represent subbands for different k± and the shaded region is the projected dispersion on k z a. 


Although we neglect anharmonicity in this paper, we hypothesize that our main conclu- 
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sions should hold even when anharmonicity is present. Phonon-phonon interaction enables 
inelastic transmission of phonons at the interface. However, the transverse symmetry selec¬ 
tion rules for k± continue to hold. Therefore, phonons crossing an abrupt or uniform interface 
have to conserve k± and are restricted to inelastic jumps within the conserving modes M c . 
On the other hand, phonons crossing a mixed interface can jump in frequency within the 
minimum of the contacts’ modes M m i n . Thus, we expect a similar relation between the 
crystal structure, which determines the ratio and the relative magnitude of the 

conductance for the abrupt, uniform and mixed interfaces. Further studies are required 
to evaluate the validity of this hypothesis. We also expect a conductance increase for all 
the systems considered in this work, since anharmonicity allows transmission of phonons 
with frequencies beyond the elastic limit. Thus as temperature and thereby anharmonicity 


increases, interfacial thermal conductance increases 


m 


26} . For some systems with large 


Debye temperatures anharmonicity can be neglected. For instance, Si/Ge interfaces present 
a thermal conductance relatively independent of temperature below 500 K, which indicates 
that phonon-phonon interactions are not dominant over that temperature range 26]. 


IV. GAIN VS. LOSS IN TRANSMISSION 

An essential part of the conductance inequality (Eq. [5j) is the transmission, which can 
be characterized from our knowledge of M and MT. For example at low frequencies, M c & 
M min and MT abr & MT uni ~ MT mix , so the transmissions are similar and they only depend 
on the acoustic mismatch between contacts. Unfortunately most of the spectrum is outside 
this low frequency regime. 

For the crystals with scalar IFCs over the mid-frequency range, MT uni > MT mix > MT abr 
as long as M c & M m , n (Fig. [3b and HJi) . Therefore, the transmission loss due to disorder for 
modes that conserve k± dominates the transmission gain from modes that do not conserve 
k±. As frequency increases, thermal energy is carried by shorter wavelength phonons and 
disorder back scattering accentuates. At some point, it becomes strong enough to reduce 
MT mix even below MT abr (Fig. 0b). 

For the DC crystal, the different polarizations available influence the transmission func¬ 
tion. For instance, mixing facilitates transmissions between TA-LA modes (Fig. [5]). This 
follows from MT mix > MT uni in regions where M min ~ M c , which implies T mix > T uni . 
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Note that this is the opposite of what we saw for the scalar crystals. Another interesting 
example shows shifting of transmission between polarizations. Around oj ~ 5.2 x 10 13 rad 
s _1 , MT mix for LA-LO decreases while MT mix for LA-TO increases. At this frequency the 
ratio M min /M c for LA-TO is larger than for LA-LO so MT mix shifts towards the more fa¬ 
vorable condition while conserving energy. In this paper we just scratch the surface of the 
importance of understanding polarization for interfacial transport, and further studies are 
required in the topic. 

To gain further insights into the transmission, we focus on the crystals with scalar IFCs. 
For uniform interfaces, Fig. [2] shows a conductance maximum when the junction mass is 


ows from a generalization of the 
271 . [28]. By Fourier transforming 


the arithmetic mean (AM) of the contact masses. This fol 
same result in ID interfaces with a single atomic junction 
the transverse coordinates, our 3D problem decouples into a sum of ID chains with IFCs 
that depend on the transverse wavevector. For each k± we assume an incident, reflected 
and transmitted wave and find their amplitudes by solving the equation of motion for the 
interfacial atom. The transmission T k± follows from the ratio of transmitted over incident 
current. In this way, MT for the uniform interface is 

MT uni = , (8) 


with 


j , _ 4r fcx r fc ^ _ 

u; 4 Am 2 + (rj^ + r£J 2 

is the broadening matrix in NEGF formalism, which reduces to a scalar function when 
dealing with a single degree of freedom per atom 22]. This quantity is related to the escape 
rate of a phonon into the contacts and is given by T k± = 2ujpv k± , with p the mass density 
and Vk ± the frequency dependent group velocity of the mode or subband defined by k±. The 
superscript in r*. x refers to the light (/) and heavy (h) contacts. Am = rrij — (mi + m ^)/2 
measures the deviation of the junction mass rrij from the AM of the contact masses. Thus, 
when rrij is the AM, each T k± , MT uni and G q uni are maxima. Note that G q uni > G q abr as long 
as mi < mj < m k , since the abrupt interface is recovered when mj = mh- Also note that at 
low frequency, w 4 Am 2 << the transmissions only depend on the acoustic mismatch 

between contacts. 

A similar generalization from its ID counterpart { 27 , 28] leads us to conclude that in an 
abrupt interface where interfacial bonding is the only variable, conductance is maximized 
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when the force constant is the harmonic mean of the contact force constants. In the same 
fashion, we can generalize other ID results to 3D interfaces 29]. 


The conductance maximum derived from Eq. [U] is not valid for tensorial IFCs (Fig. [21) ■ 
In that case, the amplitudes of the incident and transmitted waves are related through a 
matrix equation (Eq. IC16|) . Am 7 ^ 0 affects both the denominator and the numerator of 
the transmission, and therefore there is no clear trend when decreasing Am. For instance, 
Am 7 ^ 0 might abate the transmission for some polarizations but enhance the transmission 
between others. 

For mixed interfaces, we can approximate MT m i X starting from Eq. IC18I (Appendix 1C 2|) . 
the relation between incident and transmitted wave amplitudes at the interface. The heart 
of the approximation lies on Ending the inverse of the matrix (A + Zb ~ Z c ) 1 , which is a 
diagonal matrix with tiny off-diagonal elements. These small elements come from Fourier 
transforming the random mass distribution at the interface. We assume that all these 
elements are constant, since a random mass distribution contains components in the entire 
frequency spectrum. Then we estimate their value relating the known real power spectrum 
with the k space spectrum through Parseval’s theorem. Finally, we hnd the desired inverse 
using a first order Taylor expansion (( A-\-B)~ l ps A~ l — A~ 1 BA~ 1 ). With this information, 
the sum of the transmissions becomes 


MT mix = ^ T k± , k± + T k±:k ' ± , 

k L 


4TLr£, 


'-k x ,k . 


and 


Ti 


cu 4 (l — a)a(mi — m r ) 2 


k±,k' 1 


NT l k r h ki 


-Tk ± ,k ± Tk' u 


( 10 ) 


( 11 ) 


( 12 ) 


(Am) is the average over the junction masses, N is the number of atoms in the cross 
section and a is the fraction of heavy atoms at the interface. Equation [T2] suggests that 
the transmission between modes that do not conserve k±, T k± j .>^, is proportional to the 
square of the difierence between the atomic masses of the contacts, (mi — m r ) 2 , to the alloy 
scattering factor, (1 — a)a, and to some function of the acoustic properties of the contacts. 
The equation does not capture the decrease in transmission among the modes that conserve 
k± due to disorder. It also over predicts the contribution from transmissions that do not 


14 











conserve k± and does not capture their asymmetric bias as a function of junction mass 
(Fig. [ 21 ). In spite of that, it provides a sense for the expected conductance enhancement by 
mixing and insight on how to build the transmission between different modes, which is an 
important step forward towards qualitative understanding of interfacial conductance. 


V. CONCLUSION 

In this manuscript we quantify the role of crystal structure and interface morphology on 
the interface thermal conductance. We show that the crystal structure (SC: simple cubic, 
FCC: face centered cubic, or DC: diamond cubic) determines the relative magnitude of the 
minimum of the contacts’ modes M min vs. the conserving modes M c that conserve the 
component of phonon wavevector transverse to the interface k±. On the other hand, the 
interfacial morphology (abrupt, uniform: with an added homogeneous layer, or mixed: with 
atomic disorder) determines if phonons can move through M c or M m i n . 

We find that adding a unit cell monolayer of mixing to an abrupt interface enhances the 
interfacial conductance, but the degree depends on the ratio M min jM c . In particular for 
a scalar FCC crystal, where M min M c , the conductance of a mixed interface increases 
relative to that of an abrupt interface by 13%. This modest enhancement comes from a 
balance between the new accessible modes and the extra scattering created by disorder. For 
a SC crystal, the relative conductance increment from the abrupt interface to the mixed 
interface is ~ 150%. This large enhancement comes from a region where there are available 
modes but they do not overlap ( M min > M c = 0) because of a central void of modes in the 
Brillouin zone. The void appears when the upward shift in the subbands dominate their 
shrinking as k± increases, which happens due to the independence of interlayer coupling 
from k±. 

For a DC crystal, we find that the effect of mixing depends on the polarization. In par¬ 
ticular, mixing increases transmissions between TA branches but not between LA branches. 
The modes for the TA branches present a central void similar to what we saw in the SC 
crystal. This suggests that materials with modes containing central voids are prone to high 
conductance enhancement by mixing. We also find that the conductance across a mixed 
interface does not always increase relative to that at a uniform interface. In fact, it increases 
for SC and DC but decreases for FCC, and is once again correlated with the ratio M m i n /M c . 
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This suggests that the commonly invoked virtual crystal approximation alternatively over¬ 
estimates or underestimates the effect of interfacial mixing on thermal conductance. 
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Appendix A: Simulation details 


Each interface consists of two contacts joined by a layer of primitive unit cells (Fig. [[]). 
We find the interfacial thermal conductance using NEGF and assume that the crystal struc¬ 


ture, lattice constant a and 
used simplification 


25 


30 


32 


FCs are invariant throughout each system. This commonly 
provides an easy way to study thermal conductance through 


vibrationally mismatched interfaces. Moreover, the simplification is well suited for Si/Ge 
interfaces because the IFCs of these materials are very similar [1911 and therefore the differ- 

n 

ence in atomic mass is a dominant scattering mechanism |24J. The ratio between the atomic 
masses of the contacts is 3 (mi = 40 amu and = 120 amu) for all the systems but the 
diamond crystal, where we use the masses of Si and Ge. 

The IFCs for the scalar SC and FCC interfaces are built considering only nearest neighbor 
interactions described by a force constant of 45 N / m. Assuming a = 5 A, the thermal 
conductance for the abrupt interface is given by G a b r = 7.5 MW m -2 K" 1 for SC and 
G a br = 44.3 MW m" 2 K” 1 for FCC at a temperature of 300 K. Note that the value for FCC 
is ~ 6 times larger than for SC because the FCC crystal has twice the number of atoms per 
cross sectional area and its MT is ~ 3 times larger (Fig. [3}i and 0k). 

For the FCC LJ interfaces, the IFCs are extracted from the Lennard-Jones potential 
using e = 0.0503 eV, a = 3.37 A and a cut-off distance of 2.5cx. This potential generates 
interactions up to fifth-nearest atomic neighbors and corresponds to an equilibrium lattice 
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constant of a = 5.22 A. The potential is chosen to be identical to that used by English et al. 


32] to have a point of reference for benchmarking. In fact, we checked the consistency of 


our IFCs by comparing the phonon dispersions and densities of states against the reference. 
The conductance for the abrupt interface is G a br = 57.8 MW rn~ 2 K” 1 at a temperature 
of 147 K. Our non-equilibrium molecular dynamics (NEMD) simulations predict a larger 
G a br — 97.41 MW m -2 KA 1 at a temperature of 30 K due to anharmonic transmission of 
phonons beyond the cut off frequency of the heavy material. Note that very low temperature 
NEMD results, which are classical and mostly harmonic, should tend to high temperature 
NEGF results, where the Bose-Einstein distribution approaches the classical limit. 

For the DC crystal we use the IFCs from silicon extracted using Quantum Espresso, 
which is a software package for performing calculations using density functional perturba¬ 
tion theory, that has successfully predicted and matched experimental Kapitza conductance 
and thermal conductivity without any fitting parameters [33]. In this calculation, we used 
local density approximation (LDA) of Perdew and Zunger |34j with direct fitting. The cut¬ 
off energy for the planewave kinetic energy is 30 Ryd, while the k sampling is 4 x 4 x 4 
with Monkhorst-Pack method. We also considered 4 x 4 x 4 q points when calculating the 
dynamical matrix. The lattice constant for silicon is found to be 5.398 A. Our parameters 
were chosen after carefully satisfying convergence tests, and the dispersion of silicon matches 
the experimental data quite well. For simplicity in the calculations, we only consider in¬ 
teractions up to the second nearest neighbor. Our simulations predict G a br = 242.5 MW 
rn~ 2 K” 1 at 300 K. To check our code, we simulate the same interface with IFCs extracted 
from Stillinger-Weber potential and obtain G a b r = 276.6 MW rn~ 2 K _1 at 300 K, which is 
comparable to the G a b r = 280 MW m -2 K” 1 reported by Tian et. al. 19] for the same 


interface, potential and temperature. Those values are within 15% of the ones obtained 
using lattice dynamics and NEMD calculations G a br ~ 310 MW m -2 K _1 at 300 K jl6, 26]. 


All our MT calculations are done in transverse wavevector space (&q_-space) to simplify 
the 3D problem into a sum of ID independent problems. For the abrupt and uniform 
interfaces each ID chain consists of primitive unit cells. For the mixed interfaces we increase 
the size of the unit cell and randomly choose the atoms at the junction layer according to 
the desired average mass. The unit cell for SC has 36 atoms, for FCC has 32 atoms, for FCC 
LJ has 18 atoms and for diamond has 36 atoms. For the diamond crystal we also simulated 
16 and 64 atoms and did not see appreciable changes in the results. Based on this, the 
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results for FCC LJ might change less than 5% if we increase the number of atoms. The 
MT for scalar SC and FCC agree with the MT obtained using Eq. IC181 For each mixed 
interface, we report the average over more than 12 independent calculations and in Fig. [2] 
we also report the standard deviation. 

To split the contribution of MT m i X from the modes that conserve and do not conserve 
(dashed line in Fig. |3jr and 0k), we find the transmission directly from Eq. IC18I in a 
system with 40 x 40 atoms in the cross section and periodic boundary conditions. Our 
results show the average over 12 independent simulations of random distributions of atoms 
at the junction. 

To calculate propagating modes for a contact we simulate an “interface” where the leads 
and junction are the same material. In this case T = 1 because there is no interface and 
MT = M. The dispersions in the scalar SC and FCC crystals are simple enough that we 
found the propagating modes analytically by projecting the frequency iso-surface onto the 
k x , k y plane. 


Appendix B: Polarization-Resolved Transmission 


To find the transmission resolved by polarizations we start by 1) choosing a frequency 
u for which we identify all the propagating and evanescent modes of both contacts. This 
is done by solving a generalized eigenvalue problem as explained by Wang et. al. (Sec. 


2 .2.2 of 


2.3]). Then we 2) assign a polarization to each of the propagating modes. That 


is, we find the dispersion branch to which each mode belongs. This is done by moving 
in small wavevector increments from k 0 , a fixed wavevector where we know the correspon¬ 
dence between eigenvalues (frequencies), eigenvectors (polarizations) and branches, to k !, 
the wavevector of the phonon we want to label. In each step we calculate the eigenvectors 
of adjacent k grid points and project ones into the others. Then according to the maximum 
projection between eigenvectors we assign a branch to each of the eigenvalues and eigenvec¬ 
tors of the next grid point. Once we assigned a label to each propagating mode we 3) find 
the response around the interface to an incident mode from the left contact. This is done 
using the Green’s function of the system, which is the impulse response of the system, by 
exciting the system with a superposition of impulses that resemble the mode. Then we 4) 
project the part of the response at the right contact onto the modes of that contact. At 


18 







X-i Xo Xi 
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FIG. 8. System split into sites in the transport direction. 

this point we have the amplitude of the impinging and transmitted modes. Finally, we 5) 
find the current carried by each mode and the transmission between modes, which we label 
according to the labels of the modes involved. For mixed interfaces we have to unfold the 
branches of the superccll to be able to identify the polarization and label consistently with 
primitive unit cell polarizations. 

Appendix C: Transmission for Mixed Interface 

Our aim is to solve the scattering problem of a wave impinging on an interface to obtain 
Eq. [TOj [H] and [121 To this end we 1) assume incident, reflected and transmitted waves and 
find an equation relating their amplitudes. Then we 2) approximate that equation to find an 
analytical solution. Finally we 3) End the transmission from the ratio between transmitted 
and incident currents and sum them up to get MT. 

1. Equation Relating Amplitudes 

Consider a system split into sites in the transport direction (Fig. [8]) and described by the 
equation of motion 


^ MsysX ^'sysXi 


(Cl) 


with M sys the mass matrix of the system 


/• 
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Mi 


M, 


sys 


m 2 


(C2) 


M 3 
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F sys the force constant matrix of the system 




F = 

± sys 
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F on 1 F c 


off 1 


\ffi F on 2 F off3 
r+ 

off 3 


Fit, F on3 


V 


(C3) 
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and x the vector containing the displacements from equilibrium of each of the atoms of the 
system. The equation of motion for the interfacial site is given by 


u) z M 2 Xo — F on 2 X 0 + FoffiX-i + Foff 3 X 1 


(C4) 


Because of the periodicity of the contacts, plane waves of the form Xn = Xje lkjna (Bloch 
states) satisfy the equation of motion for any contact site if Xj, the polarization vector, 
satisfies 


u z MXj = 


F -l- pl P ~ ik 3 a _i_ f ,, p ik i a 


Xy 


In terms of these plane waves we assume a solution for the system of the form 


Xn 


V A k+ x k+ e lk ' na + y B k -X k - 

/ K'l K'l / -V K'l K'l 


g —na 


(C5) 


(C6) 


for n < 0, and for n > 0 




(C7) 


where + and — refer to plane waves propagating to the right or left, X and Y refer to the 
polarizations on the left and right contacts. We replace the assumed solution (Eq. IC6I and 
EZD into the equation of motion at the interface (Eq. 1( 41) . For the factor ^//iX-i, we split 
each F^^e 1 ^ 0, into Hermitian and anti-Hermitian parts. We replace the Hermitian part 


111 


favor of u 2 M\ — F on \ using Eq. IC5l and reorganize the anti-Hermitian part in matrix notation 
to get 


q w M\ — F on 1 


F off 1 X.-I - 


Xo — ZaVi+A — Z b \i_B, 


(C8) 


with Vi + and Vi_ the matrices whose columns are the polarizations X k + and X k - respectively 
and with 


f.hiWhV ',; 1 - d//iU+Arp 1+ 


1 —Itz—1 


Z A = 


(C9) 
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where the second subindex of kf x run over the possible kf. In a similar way we get that 

oj 2 M 3 — F on3 


FoffiXi — 


"Xo + ZqV 3+ C 


with 


^c = 


F am V 3 + \ 3 + Vi + ' - F^Ws+A^U" 1 


A = u 2 


m 2 - 


m x + m 3 


Fqu 2 


Fonl + F on 3 


From there we can derive the coefficients for the transmitted waves 


(C12) 


(C13) 


Equating Eq. IC6I and IC7I at n = 0 and putting Eq. IC8I and 1C 121 into Eq. IC4I we get the 
following set of equations 

Xo = + V X .B = lA+C, (C14) 

A Xo = -Z a V 1+ A - ZM-B + Z c V 3+ C , (C15) 

with 


C = V 3+ v (A + Z B - ZcY 1 (Z B - Z A ) E 1+ A 


(C16) 


2. Approximation of Amplitudes Equation 


Imagine that each site on Fig. [8] consists of a cross sectional plane of atoms in the mixed 
interface (Fig. [TJd). For the SC and FCC scalar systems, the force constants are invariant 
in the transport direction, periodic in the transverse direction and scalar between atoms. 
Thus, V = V\ + = V\- = V 3+ . V is the matrix associated with a Fourier transformation into 
the transverse fc-space, whose columns are plane waves defined by k± over the N atomic 
positions r n in a cross sectional plane 
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Yn 
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e ik ± 1 r n 
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| \ 

e ik ± 2 r n _ _ . 

1 / 


(C17) 
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Using this information we simplify the relation between the impinging and transmitted waves 


(Eq. IC16|) as 


C — A + Z 


B — 


-1 


J B 


Z a )A 


(C18) 


where the tilde means the matrix in Fourier space, i.e. Z A = V' Z A V. Because of transverse 
periodicity, all the matrices in Eq. 10181 are diagonal except 


M, 


1,3 


j\r Ym n ,r 


J(k ±j -k ±i )-r 


(C19) 


For i = j the term reduces to the average if the interfacial masses 


Mo 


= (M 2 ) = (1 - oi)mi + am h , 


(C20) 


where a is the fraction of heavy atoms at the interfacial layer. For i 7 j we are calculating 
a frequency component of a random distribution of masses, which should spam over all the 
kj_ spectrum. Thus we assume that all the off diagonal components of M 2 have the same 
magnitude. We estimate the value using Parseval’s theorem, the power spectrum in real 
space and the transformation of the interfacial mass function at k± = 0 


Mjj 


(1 — a)a 

N- 1 


I mi - m h | 


(C21) 


Plugging this simplification and Zb = — Z& into Eq. IC18l our problem reduces to solve 




k ±1 e 
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-1 


C = - 2 


6 Ck ±2 




Z Ak ±2 
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'■/ 



V '7 


A, 


(C22) 


with 


Ck ± — u 2 ( (M 2 )- Y —" ) — ( Z A k ± + Z C k ± ) 


2 / (1 — a)a 

e = uj \I — -— \mi -m h \ 


N — 1 


e is small since it is inversely proportional to \/fV — 1, so we approximate the inverse of the 
matrix using the first order of its Taylor expansion (A + B) 1 « A~ l — A~ 1 BA~ 1 with A 
being the diagonal part and B the rest. Finding the inverse and solving Eq. IC22l we get that 


C = QA 


(C23) 
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with 


—2 Z Ak± 

Qln = { 

&Z Ak , . 


if / = n 


if l^n 


(C24) 


Ck ±l (k ±rL 

Where Qi n relates the amplitude A n of the n incident mode with amplitude C) of the l 
transmitted mode. 


Find the Transmission 


Now that we know the coefficients we can calculate the transmission from mode A n to 
mode Ci by dividing the transmitted by the incident current 
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with 


Here we replace 


Cfc = OJ 2 


{m n ) - 


mi + m r 


+ i 


r lkn _|_ ^rk n 


ry • 

6Ak ± ~ ^Ck ± ~ -1- 


2 2 ’ 

which is true only for the propagating modes and therefore it works only when both of the 
modes involved in are propagating, i.e. when Tm 0. 

Then the MT per unit cell is 

rrT1 i % r j k Yrk i x—"v r r 

MT VUC = — > 
puc m 


n 1 r^n 


N ^ \Ckn 

and from there Eq. [TO] |TT] and [12] follow. 
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